The effect of nonlinearity on adiabatic evolution of light 
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We investigate the effect of nonlinearity in a system described by an adiabatically evolving Hamil- 
tonian. Experiments are conducted in a three-core waveguide structure that is adiabatically varying 
with distance, in analogy to the STIRAP process in atomic physics. In the linear regime, the system 
exhibits an adiabatic power transfer between two waveguides which are not directly coupled, with 
negligible power recorded in the intermediate coupling waveguide. In the presence of nonlinearity 
the behavior of this configuration is drastically altered and the adiabatic light passage is found to 
critically depend on the excitation power. We show how this effect is related to the destruction of 
the dark state formed in the STIRAP configuration. 



The adiabatic theorem describes one of the most pow- 
erful concepts in quantum physics[l]. It states that if the 
parameters of a quantum system evolve slowly enough in 
time, the associated initial eigenstates will be preserved, 
and there will be no exchange of energy between them. 
This well studied theorem finds wide applications in di- 
verse areas of science ranging from molecular physics to 
quantum field theory, from chemistry to nuclear physics. 
A close reexamination of the adiabatic principles led to 
the discovery of Berry's geometric phase[2| - known to 
occur ubiquitously in many processes in naturejjj- Quite 
recently, quantum adiabatic methods were suggested as 
a basis for a new class of algorithms meant to address 
NP-complete problems within the framework of quantum 
computing|4]. In addition, techniques exploiting an adi- 
abatic passage provide practical approaches in achieving 
nearly complete population transfer between two quan- 
tum states 0, IE M HI. One such example of coherent 
adiabatic excitation is stimulated Raman adiabatic pas- 
sage (STIRAP) that makes use of two appropriately pre- 
pared laser pulses in order to couple two non-degenerate 
met ast able states via an intermediate level. Remarkably 
this can be achieved without any appreciable excitation 
of the intermediate state 0, IE Ej. 

One of the underlying - and sometimes limiting- as- 
sumptions of the adiabatic theorem is the presumed in- 
trinsic linearity of the system, a condition that is often 
not met under actual experimental conditions. For exam- 
ple, nonlinearity comes into play in various adiabatically 
evolving systems such as Bose-Einstein condensates in 
slowly varying potentials or fieldsJIEEjlEEEIj] an d non- 
linear optical processes [3, EH- Of course, the question 
naturally arises on how nonlinear effects may influence 
such adiabatic transfer processes [ill EH, EE E, E3 - an 
aspect that has so far eluded experimental observation. 

In this letter we consider the influence of nonlinearity 
in systems described by an adiabatically evolving Hamil- 
tonian. Experiments are conducted in a system of cou- 
pled optical waveguides, in which the distance between 
channels changes slowly along the propagation axis. Non- 



linear optical waveguides, described by the nonlinear 
Schrodinger equation, allow one to take a simple and 
direct experimental look at the interplay between adi- 
abatic evolution and nonlinearity. In addition they pro- 
vide a direct analogy with various other quantum pro- 
cesses. These include time-dependent quantum effects in 
atomic physics, Bose-Einstein condensates in time vary- 
ing traps and time dependent quantum-well potentials 
- all described in different regimes by the same evolu- 
tion equations presented here. As an example, we use a 
three-waveguide structure that reproduces the STIRAP 
process in atomic physics[18]. In the linear regime, the 
system exhibits a complete and irreversible power trans- 
fer between two waveguides that are not directly coupled, 
via an intermediate channel. Remarkably, this intermedi- 
ate waveguide carries no significant field amplitude dur- 
ing the power exchange. In the nonlinear regime, the 
adiabatic light passage is found to critically depend on 
the excitation power levels. We show how this effect is 
related to the destruction of a dark state formed in the 
STIRAP configuration!!!]. 

Consider a system of three single-mode, evanescently 
coupled nonlinear waveguides (denoted as 1, 2 and 3, see 
Fig. la). The waveguides are identical in shape and 
have a constant width along the propagation direction, 
z. However, the distances between the waveguides vary 
along the propagation. Waveguides 1 and 3 are parallel 
to each other, while waveguide 2 is oblique; it is closer 
to waveguide 1 at z = 0, and closer to waveguide 3 at 
z = L. where L is the sample length (see Fig. la). As a 
consequence, the coupling rates between the waveguides 
vary slowly along the propagation. At z = the cou- 
pling between waveguides 1 and 2 (C12) is strong, while 
the coupling between waveguide 2 and 3 (C23) is weak. 
At the output of the system (at z = L) the situation is 
reversed, i.e. C23 > C12. The coupling between waveg- 
uides 1 and 3 is practically zero in this configuration. 

The evolution of the modal amplitudes in these three 
waveguides can be described by the following set of cou- 



2 



(a) 




(b) 



0.5 











\ / 


I T 
' \ 


isJ. ^ _ 





-1 -0.5 0.5 1 
Normalized propagation distance 

Figure 1: (color online), (a) A schematic view of the STIRAP 
sample. The relative distance between the coupled waveg- 
uides (denoted 1,2 and 3) changes slowly along the propaga- 
tion axis z, resulting in slowly changing couplings rates be- 
tween the waveguides, (b) Adiabatic light passage as calcu- 
lated from Eq. (2) for a 3-core structure with a = 66m _1 , L — 
3cm, k = 600ra -1 (see text for definitions). The intensity in 
every channel is plotted as a function of normalized distance. 



pled discrete nonlinear Schrodinger equations: 
.dE n 



dz 



" PnE n + Cn,m(z)E m + T\E n \ 2 E n = (1) 



where n = 1,2,3, E n is the wave amplitude in waveg- 
uide n, f3 n is the longitudinal wave vector (propagation 
constant) for the mode or bound state in waveguide n 
and the summation is carried out on nearest-neighbors. 
The last term in Eq.(l) accounts for the nonlinear de- 
pendence of the on-site wavevector /?, where V is associ- 
ated with the Kerr nonlinear coefficient of the waveguide 
structure. This term is important only in the nonlinear 
regime and can be neglected at low light power levels. 

In the linear limit, the description of this system by 
Eq.(l) carries a perfect analogy to the STIRAP process, 
first described in the framework of atomic physics 0, @|. 
This surprising process leads to a complete transfer of 
population between two atomic levels for which a direct 
transition is forbidden, via a third level. However, in 
the adiabatic limit the intermediate level is never popu- 
lated during the process|6]. Indeed, the equations used 
to describe the STIRAP effect in atomic physics are iden- 
tical, under the rotating wave approximation, to Eq.(l) 
in the linear limit. In this analogy z replaces time, the 
amplitude in each waveguide corresponds to the ampli- 
tude in each atomic level and the coupling between the 
waveguides plays the role of the Rabbi coupling of the 
atomic energy levels caused by resonant electromagnetic 
radiation. Identical values of the parameter/5 for coupled 
waveguides represent zero detuning of the electromag- 
netic radiation from the level spacing. A linear STIRAP 
scheme was recently suggested in an optical system us- 
ing a different analogy that required an imprint of pe- 
riodic gratings or bending of the waveguides along the 
propagation axis, to introduce coupling between dissimi- 
lar waveguides (Til l2o|. An implementation using identi- 
cal waveguides and a simple geometry similar to the one 
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Figure 2: (color online). Adiabatic passage in the STIRAP 
sample, (a) Measurement of the output light distribution 
when light is injected into waveguide 1. After the adiabatic 
sweep, the light emerges from waveguide 3. (b) The same 
experiment in a truncated sample, showing that during the 
adiabatic sweep, there is significant intensity in waveguide 2. 
(c) BPM simulations of the propagation, (d)-(f) The same 
as (a)-(c), when light is injected to waveguide 3 and emerges 
from waveguide 1. In this case, during the adiabatic sweep 
the intensity in waveguide 2 is negligibly small. 



discussed here was proposed by Paspalakis|21], and re- 
cently implemented in the linear regime by Longhi and 
coworkers [18J. 

To theoretically analyze the linear response (7 = 0) of 
the system shown in Fig. la we recall that the coupling 
coefficient between two evanescently coupled waveguides 
varies exponentially with the separation distance|22|. As 
a result, for a structure of length L, the two coupling 
constants are found to vary according to C\2(z) = k • 
exp[— a(z — L/2)] and C2z{z) = K-exp[a(z — L/2)], where 
k is the coupling strength in the middle of the structure 
(z = L/2) and a is a slow adiabatic parameter related to 
the slope of waveguide 2, that is 7 = ol/k <C 1. If at the 
input of this system, the third waveguide is excited, i.e. 
£3(0) = 1 , then by employing WKB expansion methods 
one can show that to a very good approximation the field 
in the first waveguide evolves according to: 
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In Eq.(2), A' 1 = [-4 7 2 - 2 7 2 tanh(2£ ) - 2 cosh(2t )], 
t = aL/2, t = a(z - (L/2)), -t < t < t , tanh(cV) = 
-7(2/cosh(2t )) 1/2 , and Q(t) is a phase function. E2 
and E3 are obtained by plugging Eq. (2) into Eq. (1), 
and using the conservation law \Es\ 2 — 1 — |£a| 2 — |^2| 2 - 
Fig. lb shows the evolution of the intensities I n = |i? n | 2 
in a 3-core adiabatic system with parameter values very 
close to those used in our experiments, as obtained from 
the analytical expressions of Eq. (2). The numerical 
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results are not shown here since they are very close to 
those already depicted. As clearly shown in Fig. lb, 
the power adiabatically leaves channel 3 and eventually 
populates channel 1, with very little energy remaining 
in the intermediate waveguide. This is in perfect anal- 
ogy to the STIRAP process. The first term on the right 
of Eq.(2) is primarily responsible for this adiabatic tran- 
sition whereas the second one describes the oscillatory 
component in Fig. lb. 

The waveguide triplet used in our experiment was fab- 
ricated on an AlGaAs substrate, using standard pho- 
tolithography techniques [23]. The waveguides have a 
width of 3 /im, and the sample length is L=18mm. The 
edge-to-edge distance between waveguide 1 and 2 is 2 
jam at z=0, and 7 jam at z=L, while the distance be- 
tween waveguide 1 and 3 is fixed at 12 jam. This yields 
a coupling of 2500 m _1 between waveguide 2 and 3 at 
z=L=18 mm and a coupling of 250 m _1 between waveg- 
uide 2 and 3 at z=0, while the coupling between the 
waveguides is 790 m~ l at z=L/2. A second sample with 
similar parameters was fabricated, and was truncated to 
enable observation of the amplitude in the waveguides 
before the full sweep is achieved. In the experiments pre- 
sented below, light is injected into one of the waveguides 
in the structure at z=0, propagates along the sample and 
is measured at the sample output. Nonlinearity is intro- 
duced by increasing the power of the input beam. A 
full description of the experimental setup can be found 
elsewhere [23]. 

Fig. 2 shows the result of experiments done at low 
powers. When the input beam is launched into waveg- 
uide 1 (Fig. 2a), the output light emerges from waveg- 
uide 3. However, a similar experiment done in the trun- 
cated sample (Fig. 2b), reveals that waveguide 2 carries 
a significant field amplitude during the power exchange 
between waveguide 1 and waveguide 3. This is also illus- 
trated in the BPM simulation shown in Fig. 2c. On the 
other hand, when light is initially injected into waveguide 
3, it emerges from waveguide 1 as shown in Fig. 2d, yet 
the truncated sample shows that in this case the inte- 
sity in waveguide 2 is negligible during the process (Fig. 
2e). This result is corroborated by the BPM simulation, 
as shown in Fig. 2f. Even though the coupling between 
waveguide 3 and waveguide 1 is zero, power is fully ex- 
changed between them during the adiabatic sweep, with 
only minimal excitation of waveguide 2. 

We now turn to the effect of nonlinear perturbations on 
the adiabatic passage described above. For this purpose 
we again launched light into waveguide 3, and measured 
the output light distribution as a function of the input 
beam power. The results of this experiment are presented 
in Fig. 3a. These results show that the presence of non- 
linearity reduces the efficiency of the adiabatic passage, 
even at relatively weak powers. The experimentally mea- 
sured light distribution at the output are compared to 
BPM numerical results in Fig. 3b, taking into account 
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Figure 3: (color online). The effect of nonlinearity on adi- 
abatic passage, (a) Measurements of the output light dis- 
tribution as in Fig. 2d, at different input intensities, (b) 
comparison between the experimental results (markers) and 
numerical calculations (lines, see text), (c) Numerical calcu- 
lations of the intensity distribution in the sample along the 
propagation, for an input power of 350W. 



corrections due to dispersion and cross-phase modulation 
effects [24]. The experimental and numerical results show 
good agreement in the weak nonlinear regime, while at 
higher powers the experiment deviates from the theory, 
probably due to nonlinear absorption effects. Fig. 3c 
shows an example of the calculated evolution of the in- 
tensities in waveguides 1 and 3 along the propagation in 
the nonlinear regime (power of 350W) . This figure should 
be compared with the linear dynamics in Fig lb. 

These results are compatible with previous theoretical 
predictions that considered the mean-field dynamics of a 
Bose-Einstein condensate in a time dependent triple- well 
trappjj. The authors have shown that the adiabatic pas- 
sage should break down when the magnitude of the non- 
linear parameter T exceeds that of the detuning between 
levels. In the optical analogue, detuning is introduced 
when the waveguides have different propagation parame- 
ters (3. In the configuration used here all three waveguides 
are identical, which corresponds to zero detuning, hence 
the adiabatic passage is expected to break down even for 
very weak nonlinearity. 

The STIRAP effect relyes on the existence of a dark 
eignstate of the system, a phenomenon known as Co- 
herent Population Trapping (CPT)0, It has been 
theoretically shown that dynamical level shifts induced 
by nonlinearity can affect the resonance condition that 
leads to the CPT state, hence reducing the efficiency of 
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Figure 4: (color online), (a) A schematic view of the sample 
used to probe the dark state, (b) Formation of the dark state 
in the linear regime. Light is injected to waveguide 3 and re- 
mains trapped there, despite the coupling between waveguide 
3 and waveguide 2 (see text), (c) Partial destruction of the 
dark state by nonlinearity. 

STIRAP[12]. To demonstrate this effect in our system, 
we consider the configuration presented in Fig. 4a which 
is identical to the configuration of the STIRAP sample 
at z=0, but with no variation of the couplings along the 
z direction. Waveguides 3 and 2 are weakly coupled, 
therefore light injected into waveguide 3 is expected to 
tunnel along the propagation to waveguide 2. However, 
the strong coupling between waveguide 2 and waveguide 
1 results in two new modes with propagation constants 
that are spaced symmetrically around that of the third 
guide. This leads to a sharp resonance that eliminates 
the tunneling, and therefore light that is injected into 
waveguide 3 remains trapped in that waveguide. The for- 
mation of this dark state is experimentally demonstrated 
in Fig. 4b. When nonlinearity is introduced by increas- 
ing the input power (300W), the eigenvalue of the mode 
in waveguide 3 is shifted and the resonance condition is 
no longer satisfied. As a result the dark state is destroyed 
and tunneling out of waveguide 3 is partially recovered 
(Fig. 4c). Since the STIRAP effect is based on the evo- 
lution of this dark state, this also explains the sensitivity 
of STIRAP to nonlinearity. Is it interesting to note that 
even though the level detuning due to nonlinearity can in 
principle be compensated by the sample design, the dark 
state may still be dynamically unstable[12l|. 

In summary, using coupled nonlinear optical waveg- 
uides we have investigated the effect of nonlinearity on 
an adiabatic process - an optical analogue of the STI- 
RAP process. In the nonlinear regime, we found that 
even weak nonlinearity is enough to impair the efficiency 
of STIRAP. This was explained by the destruction of the 
dark state formed in the STIRAP configuration. 

The approach presented here can be extended to more 
complex structures, implementing a variety of slowly- 
varying photonic potentials and giving rise to new non- 
linear effects. Waveguide lattices can be used to adi- 
abatically introduce changes in the dispersion relation, 
for example by opening or shifting gaps in the spectrum 
or by introducing disorder, offering a new experimental 



playground for the study of the interplay between non- 
linearity and adiabaticity. 
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